Electronic Thermal Conductivity of Multi-Gap Superconductors with Application to 

MgB 2 
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The remarkable field dependence of the electronic thermal conductivity observed in MgB2 can 
be explained as a consequence of multi-gap superconductivity. A key point is that for moderately 
clean samples, the mean free path becomes comparable to coherence length of the smaller gap over 
its entire Fermi surface. In this case, quasiparticle excitations bound in vortex cores can easily be 
delocalized causing a rapid rise in the thermal conductivity at low magnetic fields. This feature is 
in marked contrast to that for anisotropic or nodal gaps, where derealization occurs only on part 
of the Fermi surface. 

PACS numbers: 74.20.-z,74.60.Ec,74.25.Fy 



The unexpected discovery of superconductivity in 
MgB 2 with a relatively high T c = 38 KEj aroused great in- 
terest and was soon followed by experiments which estab- 
lished phonon mediated s-wave superconductivity, e.g., a 
B-isotope effecto, a coherence peak in 11 B nuclear re- 
laxation rateEl and aupftxponential dependence for tem- 
peratures T < 10 Kcn3. Theoretical studies concluded 
that the coupling of the holes in the 2p cr -bands of the 
B-planes to bond stretching modes was strong and pri- 
marily responsible for superconductivity. The electron- 
phonon coupling on the parts of the. Fpixni surface asso- 
ciated to 2p 7r -bands is much weaker&ErBQ. 

Despite its standard origin, superconductivity in MgB 2 
has several unusual properties pointing towards a more 
complex nature. One aspect is the presence of two 
gaps of different magnitude. Their ratio is esti- 
mated as r n ^p-A§iAq ~ 0.3-0.4 based on various 
experimentsaHMO£3o. Evidence for two gaps is also 
provided by the rapid rise of the specific-heat coefficient, 
j s (H), at very low magnetic fieldsoE-i Orbital depen- 
dent superconductivity has been proposed theoretically 
by several authors with the pri mar y (secondary) gap as- 
sociated with the a- (7r)-bandstall3. 

Recent studies of the inplane thermal conductivity in 
a magnetic field show an unusual field dependenceO. For 
fields both parallel and perpendicular to the c-axis, the 
electronic thermal conductivity k s (H) exhibits a steep 
increase in the low-field region, suggesting a large release 
of mobile quasiparticles in the mixed state. This con- 
trasts strongly with the behavior of conventional s-wave 
superconductors, where quasiparticles bound in the vor- 
tex cores gi ye .| V iary little contribution to n s except very 
close to iJ C 2c§li3. At first glance, a small secondary gap 
Ag in multiband models would provide enough carriers 
for transport at low fields. However, they would be non- 
mobile carriers inherent in their s-wave character. It is 
our aim here to reexamine thermal transport for multi- 
gap superconductors and show the drastic influence of 
sample purity on the characteristic behavior of k s (H). 

The measured MgB 2 samples are regarded as being in 
the moderately clean regime: experimental estimates of 
the mean free path give £ ~ 500-800 A, compared with 



the inplane coherence length £ a b ~ 120 A derived from 
H C 2, whichas determined by the primary superconduct- 
ing cr-bandEI The relevant length scale here to be com- 
pared with £ is, however, £q of the passive 7r-band. Thus, 
we may consider the quasiparticles in the er-band in the 
moderately clean regime, while those in the 7r-band with 
£o ~ £o/ r can b e m marginally clean regime. The nu- 
merical calculation based on the Bogoliubov-de Gennes 
framework shows that low-energy states in the smaller 
gap are loosely bound in vortex coresEiJ. Moreover, re- 
cent scanning tunneling spectroscopit-measurements con- 
firmed a large $ ~ 500 A in MgBjil Since the magni- 
tude of the secondary gap is small all over the Fermi 
surface, we expect a distinctively different behavior of 
k s (H) compared to single-band superconductors with an 
anisotropic gap or even gap nodes. 

In view of these circumstances, we analyze the field de- 
pendence of k s and the density of states (DOS^_JV s . For 
this purpose we introducej-Besch's solutionc3i2-fE3 for the 
quasiclassical formalisrrOEZI, which is known to be valid 
for the range of purity in question. We show that the 
quasiparticle excitations in the small gap bound in the 
vortex cores can easily be delocalized in the marginally 
clean regime, causing a rapid rise at low magnetic field. 
This field dependence is definitely stronger than that ob- 
tained for any of the single-band models. On the con- 
trary, superclean samples should exhibit a behavior very 
similar to that of conventional s-wave superconductors. 

We restrict our considerations to the case of the in- 
plane thermal current with H || z. Thus inplane impurity 
scattering is the most important for the thermal trans- 
port. The scattering—matrix between a- and 7r-bands is 
assumed to be smallo because of the different parity of 
the two orbitals under the reflection z — > —z. Thus, 
we neglect interband impurity scattering completely and 
discuss contributions from each bands independently. In 
order to calculate k s (H) and N S (H), we introduce the 
quasiclassical propagators, 



g(uj n ,k,R) 
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d£f 3 G(uj n ,k,R), (1) 



where G is the Nambu-Gorkov Green's function matrix 
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with the fermionic Matsubara frequency, u> n , the center 
of mass coordinate, R, and the relative momentum k. 
T3 is the z-component of the Pauli matrices acting on 
the particle- hole space and k = k F /\k F \ is the unit wave 
vector at the Fermi surface. They satisfy the normaliza- 
tion condition g 2 — 1 and obey the Eilenberger equations 
(H = c = ks = 1 hereafter), 



-ev F -A(R))f 3 -A(k,R) 
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-iv F -V R g = 0, (2) 



supplemented by the gap and Maxwell equations. We 
introduce the gap matrix 



over the Fermi surface. Then, we can determine the self- 
consistent (t{lu) numerically. From the linear response 
of the thermal current jhi to the temperature gr. 
— VjT, we obtain the thermal conductivity tensor) 



as 



= v F N n 



dLJ (?) sech2 iw) x 



x ( kikjRe[g h (u;)]Re[T k (u;) 



(•5) 



The comparison with the simple kinetic theory defines 
the transport lifetime, Re[rf (w)]: 



A 



A(k,R) 
-A*(k,R) 



(3) 
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and the renormalized frequency, uj n = uj n + a(uj n ), where 
a is the diagonal element of the impurity self-energy de- 
termined by the Born or the T-matrix approximation in 
this study. We neglect vertex corrections. 

Instead of solving these transport-like equations self- 
consistently, wp-adopt the Brandt-Pesch-Tewordt (BPT) 
approximatior£3. In this approximation, an Abrikosov 
solution is used for vortex lattice structures and the spa- 
tial dependence of the magnetic field is replaced by the 
external uniform field H. Only the uniform component Tj 
is kept, since the higher Fourier /^-components of g(R) 
decrease rapidly as exp(— A 2 K 2 ), A = l/y/2eH being the 
magnetic length. On the other hand, the exact spatial 
dependence of the anomalous propagators is taken into 
account including the phase variation due to the vor- 
tices. Although this theory was designed to work well 
for H < H C 2 , especially in strongly type-II superconduc- 
tors like MgB2, a detailed comparison to numerical solu- 
tions yields good agreement both for s- and d-wave su- 
perconductors over almost the whole field rangeB. This 
numerical stu 
Volovik-theor 1 



also shows that the frequently applied 
yielding 7 S cx H for an s-wave gap and 
7 S oc s/H for gaps with lines of zeros, is restricted to the 
very low-field region. This indicates the importance of 
quasiparticle transfer ^between vortices even in the rela- 
tively low-field regionB. 

By means of BPT, the solutions in eq. (|^) can be 
obtainedB formally (after analytic continuation) as 



g k (u) = 1 - iV^[2AA(k)/v F± (k)] 2 W'(u) 



-1/2 



(4) 



where u — uj[2A/v F j_(k)], W(u) — e~ u erfc(— iu) and 
uj = uj + ia. Here A(fc) denotes the spatial average of the 
gap and v F j_(k) is the component of v F perpendicular to 
the field. The real part of g^(w) is nothing but the angle- 
dependent DOS normalized by the normal-state DOS, 
N n . In order to get the closed- form solution, we use 
the Born approximation for the s-wave scattering self- 
energy, i.e., a(to) — (5£,(cl>))/2t„, where r„ is the lifetime 
in the normal state and (• • • ) represents angular average 



>7f 2AA'(fe) Re[W(u)g k (u;)} 



v F± {k) Rfi[fft(w)] 



Here scattering by the vortices appears in addition to 
quasiparticle broadening due to impurities. Note that 
eqs. (i|)-(§|) can be redu&ed to the conventional expres- 
sions in the H = limited. Moreover one finds gt = 1 
r n in the normal state. 
We concentrate on the T — > limit in this paper. The 
gap function is factorized as 



and t 



A(k)=rA f i:V / l-H/H l 



c2. 



(J) 



where r represents the smaller gap, (0 < r < 1), while 
r = 1 is used for the larger gap or the single-band case. 
The shape of the averaged gap function A(fc) is given by 
/fei e.g., /fc = 1 for an isotropic s-wave, fa = k 2 x 



7.2 



for d 



_j,2-wave and fa = 1/y 1 + ak 2 for anisotropic s- 
wavet^]. We use the square-root field dependence inferred 
from the Ginzburg-Landau theory. 

We discuss now k x s x (H) for MgB2 and the other cases 
based on this theory. For MgB2 we use for simplicity a 
spherical (cylindrical) Fermi surface for the it- (a-) band 
and the parameters n = N£/N% — 1.5, q = v F /v F = 1.5 
and r = Ag/A^ = A%/A% = 0.35. The impurity scat- 
tering rate for the a band is moderate, rj = l/2r n Ag = 
0.3. These.. pa rameters are within the range of current 
estimatesBBBH. In Fig. [l], the contribution from the 
7r-band shows a rapid rise for very low fields, while that 
from the cr-band displays rather conventional behavior. 
This rapid rise is caused by the drastic enhancement of 
the quasiparticle lifetime of the smaller gap over the en- 
tire Fermi surface as vortices are introduced. In con- 
trast, as we demonstrate for anisotropic s-wave (ani. s) 
and d x 2_ y 2 (d) in Fig. ^, the derealization of quasiparti- 
cles occurs only on parts of the Fermi surface. Here, the 
anisotropy parameter a = 15 was used. We adopted the 
unitarity limit, S = tt/2 in the T-matrix self-energy i.e., 
a = (gj,) /2r n (cos 2 S + (3^) 2 sin 2 5) for d^^-waveB. 

The sum of both bands gives k xx (H) for MgB2 in 
Fig. ^[ The overall features reproduce the experimen- 
tal data (squares)B well with the two-band model (2s) 
for 77 = 0.3. Similarly, the single-band isotropic s-wave 
model with 77 = 0.08 (s) gives a reasonable fit for Nb 
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FIG. 1: The field dependence of the inplane thermal conduc- 
tivity. For rj = 0.3, the contribution from the marginally 
clean passive 7r-band shows a rapid rise at very low field, 
while that from the active cr-band gives conventional behav- 
ior. Anisotropic s-wave and d x i _ y 2 -wave cases are given for 
comparison. 
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FIG. 2: Comparison with the experimental data of MgB2 for 
H || z (squares)LJ. The two-gap model (2s) with r\ = 0.3 (solid 
line) explains overall featpjes of the experimental darts,. The 
results for Nb (triangles)E3 are taken from Fig. 2 ofij. The 
two-gap model in the superclean limit (dashed line) shows 
a behavior similar to that of a conventional s-wave model 
(dotted line). 



(triang les)H The transport properties depend sensi- 
tively on the purity of samples, as we can see by consid- 
ering k xx (H) of the two-band model for the superclean 
regime, rj = 0.01 (dashed line). This shows a behavior 
similar to that of a conventional s-wave superconductor 
(dotted line). In the limit rj — > 0, putting uj = in 
eqs (0)-(^), we obtain the low-field expression for the 



7r-band as K x s x /K n = (it 3 / 2 /5V2)(q 2 r]/r 3 )H/H c2 . Thus, 
even in the case of small r the low-field dependence of 
k s remains small due to the factor rj in the numerator. 
In other words, in the low-field region the excited quasi- 
particles are almost localized in the vortex cores even 
in the case of the smaller gap. On the other hand, the 
slope of N S (Q) is considerably enhanced for small r as 
N s (0)/N n = {n 2 /8V2)(q/r)y/H/H c2 . It would be inter- 
esting to test this predicted change of behavior for k s (H) 
in high-quality samples. 

The thermal conductivity k s is governed by two char- 
acteristic quantities, the DOS and the transport lifetime. 
We analyze both here in order to elucidate the origin of 
the above behavior. The field dependence of the DOS is 
shown in Fig. ^, where all parameters are the same as 
those used in Fig. [|. The sharp rise of the DOS is con- 
sistent with experimcn|taLobsei'vations of 7s (-ff) in the 
polycrystalline samplesu'til. Even though there is a big 
difference between the two-band model and the single- 
gap models in k xx (H), N s (H) shows no drastic differ- 
ences apart from the presence of a residual DOS in the 
e?-wave case as H — > 0. 
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FIG. 3: The //-dependence of the DOS at w = 0. All pa- 
rameters are the same as those in Fig. ^. All curves except 
for s-wave are similar to each other (apart from the residual 
DOS in d-wave case). 



The appropriate measure of quasiparticle derealiza- 
tion is the transport lifetime in the plane. We discuss 
the lifetime of quasiparticles in the passive 7r band, where 
r = 0.35, Wp/fp = 1.5 and k _L z. The field dependence 
of Re[T s± /r n ] is shown in Fig. § for 77 = 0.01, 0.07 and 0.3. 
As expected, the transport lifetime changes drastically, 
if the marginally clean regime (77 — 0.3) is approached, 
showing a rapid rise in the low-field region. The enhance- 
ment of the quasiparticle lifetime occurs over the entire 
Fermi surface. In addition, the slope of the DOS is much 
enhanced as shown in Fig ||. These effects yield cooper- 
atively the steep rise in the thermal conductivity shown 
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FIG. 4: The r\ dependence of the inplane transport lifetime 
of the smaller gap with r = 0.35. The purity of the samples 
affects significantly on t s ± in the moderately clean regime. 



in Fig. ||. The superclean regime (r) — 0.01), in contrast, 
gives only a weak field dependence for low fields due to 
the quasiparticle localization in this case. 

Finally, we comment on the thermal conductivity for 
H gl H c i which is not covered by our theory. Experi- 
mentally, a sudden drop of k^-As observed as the mag- 
netic field barely exceeds H c ^]\. The reduction of k s is 
usually attributed to the decrease of the phonon contri- 
bution, since phonons are scattered by the quasiparticles 
in the vortex cores. Usually this kind of mechanism leads 



to a more gentle reduction of k s I13i13. For MgB2 the con- 
ditions are more complex. For T <C Ag all quasiparticle 
contributions are frozen out in the zero-field limit and 
they remain localized in the vortex cores for H ^ H c \. 
There is a stronger scattering of phonons from core states 
in multi-gap models. Since the core DOS is considerably 
larger for the 7r-band (DOS — £f/(?"Aq ) 2 ) than for the 
(T-bandE3. For T ~ Aq the quasiparticles in 7r-band are 
sufficiently excited to contribute to the zero-field thermal 
conductivity. When vortices appear, this quasiparticle 
contribution is also reduced by scattering at the vortices 
with localized quasiparticles in the cr-band. This effect 
in combination with the phonon effect leads to an even 
stronger drop of n s . These simple considerations of the 
multi-gap effect-are in good qualitative agreement with 
the experiment!!- 7 ! 

In summary, we have discussed the inplane thermal 
conductivity and the DOS in a magnetic field along the 
z-axis in the multi-gap superconductor MgB2. The rapid 
rise of k s (H) in the low field region is a special feature 
of a multi-gap superconductor in moderately clean sam- 
ples. Even in the presence of a small gap, we predict 
conventional behavior for superclean samples. This sen- 
sitivity to sample quality has to be carefully taken into 
account in the interpretation of thermal transport data 
in a multi-gap superconductor. 
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